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Abstract 

An extensive investigation is given for magnetic properties and phase transitions in one-dimensional 
Bethe ansatz integrable spin-1/2 attractive fermions with polarization by means of the dressed en- 
ergy formalism. An iteration method is presented to derive higher order corrections for the ground 
state energy, critical fields and magnetic properties. Numerical solutions of the dressed energy 
equations confirm that the analytic expressions for these physical quantities and resulting phase 
diagrams are highly accurate in the weak and strong coupling regimes, capturing the precise nature 
of magnetic effects and quantum phase transitions in one-dimensional interacting fermions with 
population imbalance. Moreover, it is shown that the universality class of linear field-dependent 
behaviour of the magnetization holds throughout the whole attractive regime. 

PACS numbers: 03.75.Ss, 05.30.Fk, 71.10.Pm 
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I. INTRODUCTION 



Bosons and fermions reveal strikingly different quantum statistical effects at low tem- 
peratures. Bosons with integer spin undergo Bose-Einstein condensation (BEC), whereas 
fermions with half-odd-integer spin are not allowed to occupy a single quantum state due 
to the Pauli exclusion principle. However, fermions with opposite spin states can pair up to 
produce Bardeen-Cooper-Schrieffer (BCS) pairs to form a Fermi superfluid. Quantum de- 
generate gases of ultracold atoms open up exciting possibilities for the experimental study 
of such subtle quantum many-body physics in low dimensions In this plat- 

form, Feshbach resonance has given rise to a rich avenue for the experimental investigation 

n 

of relevant problems, such as the crossover from BCS superfluidity to BEC |5J, fermionic 
superfluidity and phase transitions, among others 0,0, fl. Particularly, pairing and su- 
perfluidity are attracting further attention from theory and experiment due to the close 
connection to high-T c superconductivity and nuclear physics. The study of pairing signa- 
ture and fermionic superfluidity in interacting fermions has stimulated growing interest in 

nnfnn 

Fermi gases with population imbalance [9|, HO, Hll, H2j , i.e., systems with different species of 



n 



fermions [3| as well as multicomponent interacting fermions 



14 



[l5, Q. This gives rise 

to new perspectives to explore subtle quantum phases, such as a breached pairing phase 1?| 



and a nonzero momentum pairing phase of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states 
[rsl . Q and colour superfluids (J 

One-dimensional (ID) atomic gases with internal degrees of freedom also provide tunable 
)ody systems featuring novel magnetic properties and quantum phase 



interacting many 

transitions |20l . l2ll . |22| . |23| . |24J. Although the FFLO state has not been fully confirmed 
experimentally, investigations of the elusive FFLO state in the ID interacting Fermi gas 



with population imbalance are very promising 



25 



26 



38 



4Qj, 



41 



23, |2& |2a, 130, 



31 



32 



33 



34 



35, 



421 ] . Theoretical predictions for the existence of a FFLO state in the ID 
interacting Fermi gas has emerged by a variety of methods including the Bethe ansatz 
(BA) solution 26[], numerical methods 



44J. A powerful field theory approach 



28, 
35, 



30 



31 



32 



34, 



38 



42] and field theory [35, 



43|, 



44j was used to describe a FFLO state in the 



ID Fermi superliquid with population imbalance. Nevertheless, verification of the FFLO 
signature of polarized attractive fermions is still lacking via the Bethe ansatz solution. A 
recent thermodynamic Bethe ansatz (TBA) study of strongly attractive fermions shows 
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that paired and unpaired atoms form two Fermi liquids coupled to each other. The TBA 
equations indicate that spin wave fluctuations ferromagnetically couple to the unpaired 
Fermi sea. A full analysis of magnetic effects and low energy physics of spin-1/ 2 fermions with 
polarization in both the weak and strong coupling regimes, as well as a detailed discussion on 
the universality class of the magnetic behaviour in the whole attractive regime, are desirable 
in understanding such subtle paired states in ID interacting fermions with polarization. 

In this paper we provide an extended investigation of quantum phases and phase transi- 
tions for ID interacting fermions with polarization in the presence of an external field. We 
analytically and numerically solve the dressed energy equations which describe the equilib- 
rium state at zero temperature. We extend previous work on this model to derive higher 
order corrections (up to order I/I7I 3 ) for the ground state energy, magnetization, critical 
fields, chemical potentials and the external field-energy transfer relation. The phase dia- 
grams in the weak and strong coupling regimes are obtained in terms of the external field, 
density and interaction strength. In the strong coupling regime, (i) the bound pairs in 
the homogeneous system form a singlet ground state when the external field is less than 
the lower critical value H&, (ii) a normal Fermi liquid phase without pairing occurs when 
the external field is greater than the upper critical value H c2 and (iii) for an intermediate 
range H& < H < H C 2, paired and unpaired atoms coexist. However, for weak coupling, a 
BCS-like pair scattering phase occurs only when the external field H = 0, while paired and 
unpaired fermions coexist when the field is less than a critical field. Significantly, we also 
show that the universality class of linear field-dependent behaviour of the magnetization 
remains throughout the whole attractive regime. 

This paper is set out as follows. In section [Til we present the Hamiltonian and discuss 
the pairing signature for the ID fermions with population imbalance in the whole attractive 
regime. In section IIII} we present the dressed energy equations obtained from the thermo- 
dynamic Bethe ansatz equations (TBA) in the limit T — > 0. In section IIVI we present the 
magnetic properties for the model in the weak coupling regime. We solve the dressed energy 
equations in the strong coupling regime in section |VJ The explicit forms of the magnetic 
properties and the ground state energy are given in terms of the interaction strength, density 
and external field. In section IVI| we present the full phase diagrams for the whole attractive 
regime. Section IVHI is devoted to concluding remarks and a brief discussion. 
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II. THE MODEL 



The Hamiltonian 



47 



481 ] we consider 

L 



n 



£ 

-It/ 



(j)](x 







£ 2 d 2 



4>Ax)dx 



2m dx 2 / 
j (x) (fit (x) 0| (x) 0, (x) <ix 



(1) 



describes AT 5-interacting spin-| fermions of mass m constrained by periodic boundary con- 
ditions to a line of length L and subject to an external magnetic field H. In this formulation 
the field operators <pi and 0-f describe fermionic atoms in the respective states | |) and 
| I). The 5-type interaction between fermions with opposite hyperfine states preserves the 
spin states such that the Zeeman term in the Hamiltonian is a conserved quantity. For 
convenience, we use units of h = 2m = 1 and define c = mgm / h 2 and a dimensionless inter- 
action strength 7 = c/n for the physical analysis, where n = N/L is the linear density. The 
inter-component interaction can be tuned from strongly attractive (giD — > —00) to strongly 
repulsive (giD — > +00) via Feshbach resonance and optical confinement. The interaction is 
attractive for giv < and repulsive for g\r> > 0. 

The model ([!]) was solved by Yang [47| and Gaudin [48J in the 1960's and has receiver 



renewed interest in connection with ultracold atomic gases 
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39, 



40 



41 
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22 



23|, 



24 



25|, 



26 



27, 



42i | . The energy eigenspectrum is given by 



E 



7^ J^jli ^1 > where the Bethe ansatz (BA) wave numbers {k{\ and the rapidities {A a } 
for the internal spin degrees of freedom satisfy the BA equations 

A £ + ic/2 



N 

n 



exp(i/CjL) 

A a - kt + i c/2 
A a — ke — i c/2 



m , 

TT — 

l = l^-A £ - ic /2 



M 

n 



A a - A/3 + ic 
A n - A/3 - i c 



(2) 



Here j — 1, . . . , N and a = 1, . . . , M, with M the number of spin-down fermions. 

The solutions to the BA equations ([2]), as depicted in Figure HJ provide a clear pairing 

the model. The BA root distributions in 



the complex plane were studied recently 








the 


23 




29] 



291 ] . The ground state for ID interacting 
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fermions with repulsive interaction has antiferromagnetic ordering 



23 



49 



50|L Rather 



subtle magnetism for the model with repulsive interaction was recently studied 



51 



52j. For 



attractive interaction, fermions with different spin states can form BCS pairs with nonzero 
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31 



32 



34 



35|]. 



centre-of-mass momenta which might feature FFLO states 

In the weakly attractive regime, the weakly bound Cooper pairs are not stable due to 
thermal and spin wave fluctuations. The unpaired fermions sit on two outer wings in the 
quasimomentum space 0, Q due to the Fermi pressure (see Figured]). The ground state 
can only have one pair of fermions with opposite spins having a particular quasimomentum 
k. The paired fermions occupy the central area in the quasimomentum k space. Indeed we 
find from the BA equations ([2]) that in the weak coupling limit, i.e. L\c\ <C 1, the imaginary 



part of the quasimomenta for a BCS pair is proportional to y/\c\/L. However, for strongly 
attractive interaction, i.e. L\c\ 3> 1, the BCS pair has imaginary part ±i||c|. In this regime, 
the lowest spin excitation has an energy gap which is proportional to c 2 . For the cross-over 
regime, i.e. L\c\ ~ 1, the imaginary part ±iy is asymptotically determined by the condition 
ytanh(|yL) ~ ~|c|. For this cross-over regime, the spin gap might be exponentially small. 
However, it is hard to analytically determine this small energy gap from the BA equations 
([2]). Nevertheless, for the weak coupling regime L\c\ <C 1, the bound state has a small 
binding energy eb = f?n\^\jmL which has the same order of 7 as the interacting energies 
of pair-pair and pair-unpaired fermions. In this limit, the real parts of the quasimomenta 
satisfy the Gaudin model-like BA equations 23j, |36( which describe BCS pair-pair and pair- 



unpaired fermion scattering. They form a gapless superconducting phase. Using the above 
BA root configuration, the ground state energy per unit length is given by 231 ] 



h 2 n 3 
2m 



E 



2 v ; 



7i- 

12 



7T 
1 



P< 



(3) 



in terms of the polarization P = (N — 2M)/N. The first term in Eq. ([3]) includes the 
collective interaction energy (pair-pair and pair-unpaired fermion scattering energy) and 
the binding energy (internal energy). We see clearly that for large polarization (P « 1) 
the small portion of spin-down fermions are likely to experience a mean-field formed by the 
spin-up medium. This is consistent with the observation of Fermi polarons in an attractive 
Fermi liquid of ultracold atoms 37]. 

On the other hand, when the attractive interaction strength is increased, i.e. L\c\ 3> 1, 
the bound pairs gradually form hard-core bosons while the unpaired fermions can penetrate 
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into the central region in the quasimomentum space (see Figured]). The main reason for the 
unpaired fermions and BCS pairs having overlapping Fermi seas is that in ID the paired and 
unpaired fermions have different fractional statistical signatures such that they are allowed 
to pass into each other in the quasimomentum space. In the thermodynamic limit, i.e., 
L — > oo, N — > oo with N/L finite, the binding energy of a pair is eb = h 2 n 2r y 2 / '(4m). 
The dimensionless interaction strength 7 = c/n is inversely proportional to the density 
n. This signature leads to different phase segments in ID trapped fermions 25|, |26|, |27] 
than the phase separations in 3D trapped interacting fermions, where the Fermi gas has 
been separated into a uniformly paired inner core surrounded by a shell with the excess of 

nnn 

unpaired atoms [a, 13, 13] ■ 

From the ground state energy for the model with strong attraction and arbitrary polar- 
ization 23j, we find the finite-size corrections to the energy in the thermodynamic limit to 



be given by 

E(L, N) - LE™ « + O, (4) 

where the central charge C — 1 and the group velocities for bound pairs Vb and unpaired 
fermions v u are 

up(l - P) / (1 - P) 4P 



171 171 



Here the Fermi velocity is vp = him/m. In the above equation E^ 3 is the ground state 
energy in the thermodynamic limit 



2m y A 3 

The nature of the finite-size corrections indicate that two Fermi liquids couple to each 
other and have different statistical signatures. The low energy physics is dominated by 
the charge density fluctuations. The spin wave fluctuations are frozen out. In order to 
understand the pairing signature and the subtle FFLO states in ID, one should investigate 
density distributions, pairing correlations and thermodynamics, which we do here through 
the thermodynamic Bethe ansatz formalism. In particular, we shall focus on magnetic 
properties and quantum phase transitions for the whole attractive regime. 
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FIG. 1: Bethe ansatz root configurations for pairing and depairing in quasimomentum space. For 
weakly attractive interaction, unpaired roots sit in the out-wings due to strong Fermi pressure. 
For strongly attractive interaction, unpaired roots can penetrate into the central region, occupied 
by the bound pairs. However, for weakly repulsive interaction the roots with up- and down- spins 
separate gradually. In the strongly repulsive regime, the model forms an effective Heisenberg spin 
chain with the antiferromagnetic coupling constant J ~ —4Ep/j [521], where Ep = n 2 ir 2 /3 is the 
Fermi energy. 



III. THERMODYNAMIC BETHE ANSATZ 



The thermodynamic Bethe ansatz (TBA) provides a powerful and elegant way to study 
the thermal properties of ID integrable systems. It also provides a convenient formalism to 
analyze quantum pha se transitions and magnetic effects in the presence of external fields at 



zero temperature 



45 



53 



54 



551 ]. In the thermodynamic limit, the grand partition function 
is Z = tr(e~ r<: / r ) = e~ G ^ T , in terms of the Gibbs free energy G = E — HM Z — — TS 



and the magnetization H, the chemical potential /i and the entropy S 45 



53 



54 



55j. The 



TBA equations for the attractive regime are much more subtle and involved compared to 
those for the repulsive regime. In general the equilibrium states satisfy the condition of 
minimizing the Gibbs free energy G with respect to the particle and hole densities for the 
charge and spin degrees of freedom that generates the TBA equations (details are given in 



53l . l55l . l56j]). At zero temperature, the ground state properties are determined in terms of 



7 



the dressed energies for the paired e b and unpaired fermions e u and the function 



by 



( \ - 1 m '°' (7\ 

Ml, "^(mc/2)Hx 2 {) 



e b (A) = 2(A 2 -/I- -c 2 ) - / a 2 (A - A')e b (A')tZA' 

4 J-B 

ai {A-k)e u -{k)dk, 
1 /" B 

e u (A;) = (k 2 -fi--H)- / ai (A;-A)e b (A)dA (8) 



which are the dressed energy equations [221, |45|, |53[] obtained from the TBA equations in 



the limit T —>■ 0. The superscripts ± denote the positive and negative parts of the dressed 
energies, with the negative (positive) part corresponding to occupied (unoccupied) states. 
The integration boundaries B and Q characterize the Fermi surfaces of the bound pairs and 
unpaired fermions, respectively 

The Gibbs free energy per unit length at zero temperature is given by 

1 r B 1 r Q 

G(n,H) = - e b (A)dA + — / e u -{k)dk. (9) 

7T J-B 27T J_q 

The TBA equations provide a clear configuration for band fillings with respect to the external 
field H and chemical potential fi. The polarization P varies with respect to the external 
magnetic field. From the Gibbs free energy we have the relations 

-dG(/2,H)/dfi = n, -dG(ji,H)/dH = nP/2. (10) 

IV. MAGNETIC PROPERTIES IN THE WEAK COUPLING REGIME 

For weak coupling |c| — > caution needs to be taken in the thermodynamic limit. On 
solving the discrete BA equations ([2]) in the regime L\c\ 1 the ima gina ry part of the 
BCS-like pairs tends to v/fcf/ L 23]. However, the TBA equations 221, [53 ] usually follow 
from the root patterns in the thermodynamic limit, i.e. L, N — > oo with N/L is fixed. Under 
this limit, we naturally have the BA root patterns kj = Aj ± i\\c\ with j = 1, . . . , M for 
the charge degree and the string patterns with equally spaced imaginary distribution for 
spin rapidity A™ ■ = A^ + i|(n + 1 — 2j)c, with j = 1, . . . , n. Here the number of strings 



S 



a = 1, . . . , N n . A™ is the position of the centre for the length-n string on the real axis in A- 
space. Therefore, in the weak coupling limit, i.e., |c| — > 0, the integral BA equations and the 
TBA equations do not properly described the true solutions to the discrete BA equations (j2J) 
unless under the thermodynamic limit conditions. Nevertheless, the discrepancy is minimal, 
i.e., it is 0(7 2 ). 

The BA equations ([2]) in principle give complete states of the model. However, at finite 
temperatures, the true physical states become degenerate. The dressed energies in the 
TBA equations jSJ) characterize excitation energies above the Fermi surfaces of the bound 
pairs and unpaired fermions. All physical quantities, for example, free energy, pressure 
and magnetic properties can be obtained from the TBA equations withought deriving the 
spectral properties of low- lying excitations. In the weak coupling limit, the interaction 
energy is proportional to |c| which is much less than the kinetic energy. Therefore, in this 
regime, the exact ground state energy with leading term of order |c| is precise enough to 
capture the nature of phase transitions and magnetic ordering. 

From the ground state energy ([3]) we have the relation between the external field and 
magnetization 

#~ \2n 2 m z + 4h\m z ] (11) 

2m 

where m z = M z /n and the magnetization is defined by M z = nP/2. A linear field-dependent 
behaviour of the magnetization is observed. Figure [2] shows the magnetization vs the field H 
for different interaction values |7|. We observe that the analytic results plotted from (fTTl) are 
in excellent agreement with the numerical curves evaluated directly from the dressed energy 
equations (jSJ). We also find that a fully paired ground state only occurs in the absence of 
the external field. However, for H > H c where 

# c = nV + 2| 7 |] (12) 

the fully polarized phase occurs. Paired and unpaired fermions coexist in the intermediate 
range < H < H c . The phase diagram for weak coupling is illustrated in Figure |3j 

V. SOLUTIONS TO THE DRESSED ENERGY EQUATIONS 

In this section we solve the dressed energy equations jSJ) analytically in the strong coupling 
regime to obtain explicit forms for the critical fields and magnetic properties in terms of 
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FIG. 2: Magnetization M z /n vs external field H in units 2m = H = 1 with weak coupling |c| = 1 
for different densities n. The dashed lines are plotted from the analytic result (jlip . Excellent 
agreement between the analytic result and numerical solution of the integral equations ([8]) (solid 
lines) is seen. 

the interaction strength 7. We present a systematic way to obtain these physical properties 
up to order which gives a very precise phase diagram for finite strong interaction. Here 
we note that Iida and Wadati [29J have presented a different method to solve the dressed 
energy equations. We have solved the dressed energy equations (JH1) numerically in the 
whole attractive regime to compare with the analytic results. Excellent agreement between 
numerical and analytical results is found. 



First consider the ground state P = 0. Following the method developed in [27| where 
the dressed energy equations ([8]) are asymptotically expanded in terms of l/\c\, the ground 
state dressed energy equation for P = is given by 

£ >(A). 2 (^- P )-M£ _^l_ dA - (13) 
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mixed phase 




lcl = l 



FIG. 3: Phase diagram for weak coupling value |c| = 1. The dashed line plotted from the analytic 
result (|12p is in excellent agreement with the coloured phases which are obtained from numerical 
solutions of the dressed energy equations flSJ). 



with // = )[/ + %-. For convenience, we introduce the notation 



V 



7T 
1 

2^ 



B 



B 



\A)dA 



e u {k)dk 



(14) 
(15) 



for the pressure of bound pairs and un-paired fermions. Since the Fermi point B is finite, we 
can take an expansion with respect to A' in the integral in Eq. (|T3l) . By a straightforward 
calculation, the pressure p b is found to be 



128 



6 4 fB 3 \ 2Bp b 

7T \ 6 } 7T C 407T^ CP 3 



(16) 



We obtained this equation by iteration in terms of p b and /2. In such a way, the accuracy of 
physical quantities can be controlled to order l/\c\. This provides a systematic way to obtain 
accurate results from dressed energy equations. It is free from restriction on the integration 
boundaries B and Q. Furthermore, from Eq. (1161) and the condition e b (±B) = we find 



„9 V » , ,5 



(17) 
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Finally, using the above equations and the relation dp b / dfi = n, the pressure per unit length 
follows as 

b h 2 n 3 n 2 n 3 A3 3 1 . tt 2 . 

P " ^^r{ 1 + 2W\ + W^ + w i5 -^ )] (18) 



and the energy per unit length is 

h 2 n 3 f : 2 7r 2 



1 3 1 / 7T 



1+ | 7 | + 4| 7 | 2 + 2| 7 | 3 y 1 15 



■2 



(19) 



2m [ 4 48 

The dressed energy equations (jHJ) can also be solved analytically for < P < 1. Following 



27| . we define jl — H/2. We notice that the Fermi points Q and I? are still finite in the 
presence of an external field H. Similar to the case P = 0, using the conditions e b (±B) = 
and e u (±Q) = 0, we obtain the relations 

, 4 / B 3 \ 2Bp b 8Bp u 
V ~ — -//£ n -r- 

H \ 6 J 7T |c| 7T |c| 

128^ 64(#g)j 64(#g)j 
45vr 2 |c| 3 9tt 2 |c| 3 15tt 2 |c| 3 ' 1 ' 

Q (Q 2 ~\ 2Qp" 

7T \ 6 J 7T | C| 

, 64fr)§(fli | 64(/2)^(/2)l (21) 



and 



9tt 2 |c| 3 15tt 2 |c| 



„, / 2p" 4ui 



2|c| |c| 37r|c| 



2 



Kiflifi 4(/i)2 + 16(/2) 



37T | c| 3 157r|c| 3 
^ 2 _ 2p b 64(/Z)i/2 64(/2)l 



After eliminating I? and Q, we have 



3 

24(/2)f + 16(/t)f + 80(/i)i/i^ 2 



P 3vr I ^ 2|c| + 15vr|c| 
160(/i) 3 + 640( / u)i(/2)l 



457T 2 |c| 



(23) 



^ 2 / 2p> 64(^)1/2 64(#\ f 128(/2)f(/l)l 
P ~ 3tt I ^ |c| 3tt|c| 3 15tt|c| 3 / 9tt|c| 3 ' 1 ' 
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Obviously, the pressures p b and p u are functions of fl, fl and the interaction strength c, 
i.e., p b = p b (p,,p,, |c|) and p u = p u (fl,fi, |c|). Furthermore, taking into account the relations 
^ + ^ = P/2 and + = n, after a tedious calculation we find the effective chemical 
potentials for the pairs // = /i + €b/2 and for the unpaired fermions /i" = jl = ji + i//2. 
Explicitly, 

« ~ hW S P 2 , (1-P)(49P 2 -2P+1) 

r" ~ n I 



2m t 12 l7l 
(1 -P) 2 (93P 2 + 2P+ 1) (l-P)r 9 A A 

+- ^hri " " Wur 144177 p ~ 7950P 

87^ 240|7| J L 

-324vr 2 P 3 + 15720P 3 - 7620P 2 + 166tt 2 P 2 - 120P - 4vr 2 P - 30 + tt 2 ] } , (25) 

b frW f(l-P) 2 (3P + 1)(6P 2 -3P+1) 
V ~ ^ \ 77 + 



2m L 16 12 |7 
1 - P)(5 + 17P - P 2 + 491P 3 ^ 



+ — [15(1 + 2P 2 ) + 7470P 3 + 10vr 2 P 2 



64 7 2 ' 240 1 7 1 

-180tt 2 P 3 + 335tt 2 P 4 - 420tt 2 P 5 - 15405P 4 - vr 2 + 75P + 7815P 5 ] } . (26) 

These results give rise to a full characterization of two Fermi surfaces. The total chemical 
potential can be determined from either /i 6 = /i + £(,/2 or from //" = // + H/2. The chemical 
potentials for the fermions with spin-up and spin-down states are determined by /if = fi+H/2 
and /if = fi — H/2. The energy for the model with arbitrary population imbalances can be 
obtained from E/L = \m — G + HP/2, with result 

E h 2 n 3 n 2 f (1-P)7 2 tt 2 (1 - 3P + 3P 2 + 15P 3 ) 



L 2m { 4 48 

vr 2 (l - P)(l + P - 5P 2 + 67P 3 ) tt 2 (1 - P) 2 (l + 5P + 3P 2 + 247P 3 ) 



+ 



48|7| 647 2 

_ p) 

[-15 + 31125P 4 + 1861tt 2 P 5 - 15765P 5 - 659tt 2 P 4 
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1440|7| 

+346tt 2 P 3 - 14vr 2 P 2 + tt 2 P + tt 2 - 105P - 150P 2 - 15090P 3 ] } . (27) 
This result provides higher order corrections in terms of the interaction strength y-r compared 



to previous studies 
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29] 



VI. QUANTUM PHASE TRANSITIONS 

In section HVl we examined magnetic effects and phase transitions for spin-1/2 weakly 
attractive fermions with polarization. As the attractive interaction strength I7I increases, the 
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FIG. 4: Schematic dressed energy configuration for the gapless phase in the vicinity of H&. 

bound pairs become stable and form a singlet ground state. The ground state configuration 
is characterized by an empty unpaired Fermi sea, whereas the Fermi sea of the bound pairs 
is filled up to the Fermi surface. The first critical field value H c \ diminishes the gap, thus 
the excitations are gapless. This critical field indicates a phase transition from a singlet 
ground state into a gapless phase where two Fermi liquids of paired and unpaired fermions 
couple to each other. These configurations are depicted in Figure HI 

Analysis of the dressed energy equations (JSj) reveals that a fully paired phase with mag- 
netization M z = is stable when the field H < if cl , where 

h 2 n 2 



H c i 



2m 



(28) 



2 8 V~ 4| 7 | 2 |7| 3 / 

In the vicinity of the critical field H&, the system exhibits a linear field-dependent magne- 
tization 

trz 2(#-f/ cl )/ 2 11 81-tt 2 \ 



with a finite susceptibility 



2 / 2 11 81-7T 2 \ 
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FIG. 5: Magnetization M z jn vs the external field H/e^ for c = —10 in the units 2m = h = 1 for 
different densities n. The dashed lines are plotted from the analytic result (|3ip . The solid curves 
are obtained from numerical solutions of the dressed energy equations ([8]). Excellent agreement 
is seen between the analytic and numerical results. The inset shows similar comparison between 
analytic and numerical results for the susceptibility vs external field H/e^. 

This universality class of linear field-dependent magnetization behaviour is also found 
for the multicomponent Fermi gases with attractive interaction [5^]. However, it differs 
subtly from the case of a Fermi-Bose mixture due to the different statistical signature of a 

n 

boson and a bound pair of fermions with opposite spin states [52]. For the model under 
consideration here the magnetic properties in this gapless phase can be exactly described by 
the external field-magnetization relation 

l H = 1 -e b + ^-^ b (31) 

where fj, u and /j, b are given by (|25p and (I26p with P = 2M z /n = 2m z . This relation reveals an 
important energy transfer relation among the binding energy, the variation of Fermi surfaces 
and the external field. This relation might provide evidence for the pairing signature in a 
ID imbalanced Fermi gas with attractive interaction, i.e., pairs with nonzero centre-of- 
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A 

dressed energy 




FIG. 6: Schematic dressed energy configuration for the gapless phase in the vicinity of H, 



c2- 



mass momenta. The lower critical field is reminiscent of the Meissner effect, whereas the 
upper critical field determined by ( I3TI) is reminiscent of a quantum phase transition from 
superconducting to normal states [58|. Figure [5] shows the magnetization vs external field 
for different values of the interaction strength 7. Numerical solution of the dressed energy 
equations (jHJ) shows that the analytic results are highly accurate in the strong and finitely 
strong coupling regimes. 

A similar configuration occurs for the external field exceeding the upper critical field 
given by 



h 2 n 2 
2m 



H, 



c2 



7 _i_ 9^2 
h In 

2 



4 16tt 2 , 

^aR+iwjJ (32) 

where a phase transition from the mixed phase into the normal Fermi liquid phase occurs. 
Figure [6] shows this configuration in the dressed energy language. From the relation fl3Tl) . 
we obtain the linear field-dependent magnetization as 



M z 



n 



1 - 



H c2 — H 



4n 2 7r 2 



1 + 



4 



171 7 



12 

~~2 



16(tt 2 - 6) 



(33) 



with a finite susceptibility 

X 



4 12 
1 + ^ + ^7 

7 7 



16(tt 2 - 6) 

oi |3 

3 l7l 



(34) 
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A typical phase diagram in the n — H plane for finite strong interaction is shown in 
Figure [7J Smooth magnetization curves at the critical fields H c \ and H C 2 indicate second 
order phase transitions. Very good agreement is observed between the curves obtained from 
the numerical solution of the dressed energy equations and the analytical predictions (|28p 
and (132]) for the critical fields. 




n 

FIG. 7: Phase diagram for finite strong interaction with |c| = 10. The dashed lines are plotted 
from Eqs. (|28p and (|32p. The coloured phases are obtained by numerical solution of the dressed 
energy equations (jSJ). The numerical phase transition boundaries coincide well with the analytic 
results pS} and (pD- 

VII. CONCLUSION 

In summary, we have studied magnetic properties and quantum phase transitions for the 
ID Bethe ansatz integrable model of spin- 1/2 attractive fermions. Previous work on this 
model has been extended to derive higher order corrections for the ground state energy, 
pressure, chemical potentials, magnetization, susceptibility and critical fields in terms of the 
external magnetic field, density and interaction strength. The range and applicability of 
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the analytic results have been compared favourably with numerical solutions of the dressed 
energy equations. The universality class of linear field-dependent behaviour of the phase 
transitions in the vicinity of the critical field values has been predicted for the whole attrac- 
tive re gim e. This universal behaviour is consistent with the prediction for the ID Hubbard 
model 59|. However, it appears not to support the argument 45|, |46j for a van Hove- type 
singularity of quantum phase transition for ID attractive fermions. Finite temperature 
properties of ID interacting fermions will be considered elsewhere. 

We further confirm that ID strongly attractive fermions with population imbalance ex- 
hibit three quantum phases, subject to the value of the external field H 2^: (i) for H < H c \ 
bound pairs form a singlet ground state, (ii) for H > H c2 a completely ferromagnetic phase 
without pairing occurs, and (iii) for the intermediate range H c i < H < H c2 paired and 
unpaired atoms coexist. The typical phase diagram is as depicted in Figure However, 
for weak coupling, the BCS-like pairing is unstable. Two quantum phases emerge when the 
external field is applied: (i) a fully polarized phase for H > H c , and (ii) a coexisting phase 
of paired and unpaired fermions for < H < H c . The phase diagram for weak coupling is 
illustrated in Figure [3j We have shown that the mixed phase in ID interacting fermions with 
polarization can be effectively described by two coupled Fermi liquids. Our exact phase di- 
agrams for the weak and strong coupling regimes also provide a space segment signature for 
an harmonically trapped Fermi gas in ID geometry. These quantum phases and magnetic 
oroperties may also possibly be observed in experiments with ultracold fermionic atoms 



60,1611. 
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